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Abstract 



m 

' I The Exceptional Super symmetric (SUSY) Standard Model (EgSSM) predicts three 

O , families of Higgs doublets plus three Higgs singlets, where one family develops vac- 

uum expectation values (VEVs), while the remaining two which do not are called 
CN ' Inert. The model can account for the dark matter relic abundance if the two light- 

' est Inert neutralinos, identified as the (next-to) lightest SUSY particles ((N)LSPs), 

have masses close to half the Z mass. In this case we find that the usual SM-like 
i Higgs boson decays more than 95% of the time into either LSPs or NLSPs. The lat- 

. ter case produces a final state containing two leptons l^l^ with an invariant mass 

■ less than or about 10 GeV. We illustrate this scenario with a set of benchmark 

points satisfying phenomenological constraints and the WMAP dark matter relic 
^ [ abundance. This scenario also predicts other light Inert chargino and neutralino 

^ ', states below 200 GeV, and large LSP direct detection cross-sections close to current 

limits and observable soon at XENONIOO. 
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1. Introduction 



The discovery of the Higgs boson, the last missing piece in the Standard Model (SM) of 
particle physics, is one of the main goals of upcoming accelerators. It is expected that the 
Higgs particle will be detected at the Large Hadron Collider (LHC) in the near future. 
The strategy for Higgs searches depends on the decay branching fractions of the Higgs 
boson to different channels. Physics beyond the Standard Model may affect the Higgs 
decay rates to SM particles and give rise to new channels of Higgs decays requiring a 
drastic change in the strategy for Higgs boson searches (for recent reviews of nonstandard 
Higgs boson decays see |T]). In particular, there exist several extensions of the Standard 
Model in which the Higgs boson can decay with a substantial branching fraction into 
particles which can not be directly detected. Invisible Higgs decay modes may occur 
in models with an enlarged symmetry breaking sector (Majoron models, SM with extra 
singlet scalar fields etc.) [2]- [5], in "hidden valley" models [6], in the SM with a fourth 
generation of fermions [Z]-[B], in the models with compact and large extra dimensions [S], 
[S]-[in], in the littlest Higgs model with T-parity mj etc. 0. 

Another example wherein invisible decay modes can occur is supersymmetry (SUSY), 
with the lightest Higgs boson decaying into the lightest SUSY particle (LSP). R-parity 
conservation ensures the stability of the LSP so that the LSP can play the role of cold 
dark matter (CDM) [13]. In most scenarios the LSP is the lightest neutralino, a linear 
combination of neutral electroweak (EW) gauginos and Higgsinos. In some regions of 
the parameter space the lightest Higgs boson in the Minimal Supersymmetric Standard 
Model (MSSM) decays into the lightest neutralino with a relatively large branching ratio, 
therefore giving rise to invisible final states if R-parity is conserved . 



LEP and Tevatron data still allow the neutralino LSP to be sufficiently light that 
the decays of the lightest Higgs into these neutralinos is kinematically allowed. Light 
neutralinos can annihilate efficiently through a Z-pole resulting in a reasonable density of 
dark matter. Moreover the Cryogenic Dark Matter Search (CDMS) experiment recently 
reported the observation of 2 events possibly due to dark matter scattering with an ex- 
pected background of about 0.8 events [15]. CDMS events suggest that the mass of the 
dark matter particles are around 40 — 80 GeV while their spin independent elastic cross 
section is asi ~ few x lO"'^'* cm^. If recent results of the CDMS experiment get confirmed 
then scenarios with invisible decays of the Higgs boson will become rather plausible. 

Certainly the presence of invisible decays modifies considerably Higgs boson searches, 
making Higgs discovery much more difficult. If the Higgs is mainly invisible, then the 

^In the context of the so-called "nightmare scenario" , in which the LHC produces nothing beyond the 
SM, Higgs decays into and interactions with dark matter particles were discussed in [T^j. 
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visible branching ratios will be dramatically reduced, preventing detection in the much 
studied channels at the LHC and the Tevatron. In the case where invisible Higgs decays 
dominate it is impossible to fully reconstruct a resonance and it is very challenging to 
identify it at the collider experiments, i.e. quantum numbers remain unknown. At e^e~ 
colliders, the problems related to the observation of the invisible Higgs are less severe 
[3], [16] since it can be tagged through the recoiling Z. Presently, the LEP II collaborations 
exclude invisible Higgs masses up to 114.4 GeV [izj^ . 




On the other hand, Higgs searches at hadron colliders are more difficult in the pres- 
ence of such invisible decays. Previous studies have analysed ZH and WH associated 
production [1], [I8]-[I9] as well as tiH production [20j and ttVV {bhVV) production [21] 
as promising channels. The possibility of observing an "invisible" Higgs boson in central 
exclusive diffractive production at the LHC was studied in |8j. Another proposal is to 
observe such an invisible Higgs in inelastic events with large missing transverse energy 
and two high Et jets. In this case the Higgs boson is produced by VV fusion and has large 
transverse momentum resulting in a signal with two quark jets with distinctive kinematic 
distributions as compared to Zjj and Wj] backgrounds [ID] , [IS] , [22] • 

Consideration of the possibility that the dominant Higgs decays will be invisible would 
lead to drastic changes in the strategy of Higgs boson searches. Therefore, it is a rather 
interesting subject of investigation as to the nature and extent of invisibility acquired by 
Higgs, and how it can be related to specific aspects of the models concerned, especially 
well motivated SUSY extensions of the SM. In this article, we consider the exotic decays of 
the lightest Higgs boson and associated novel collider signatures within the Exceptional 
Supersymmetric Standard Model (EgSSM) [23]- [21]. This inspired SUSY model is 
based on the low-energy Standard Model gauge group together with an extra U{1)]^ 
gauge symmetry under which right-handed neutrinos have zero charge. In the EgSSM the 
problem is solved in a similar way as in the Next-to-Minimal Supersymmetric Standard 
Model (NMSSM), but without the accompanying problems of singlet tadpoles or domain 
walls. Because right-handed neutrinos do not participate in the gauge interactions in this 
model they can be superheavy, shedding light on the origin of the mass hierarchy in the 
lepton sector and providing a mechanism for the generation of the baryon asymmetry in 
the Universe via leptogenesis |25j . 

The particle spectrum of the EgSSM contains exotic matter. In particular, it involves 
three SM singlet superfields that have non-zero U{1)n charges. One of these singlets 
acquires a non-zero vacuum expectation value (VEV), breaking U{1)n symmetry and 

^Similar limits could apply to the case where the Higgs decays some fraction of the time into soft 
lepton pairs plus missing energy, as will be the case for some of the novel Higgs decays considered in this 
paper. 
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inducing the effective term and the masses of the the exotic fermions. The masses of 
the fermion components of the two other singlet superfields (Inert neutrahnos) are also 
related to the VEVs of the Higgs doublets. Because the Yukawa couplings that determine 
the strength of these interactions are constrained by the requirement of the validity of 
perturbation theory up to the Grand Unification scale the masses of the corresponding 
Inert neutralinos are expected to be lighter than 60 — 65 GeV. As a result the lightest 
Inert neutralino tends to be the lightest SUSY particle in the spectrum. Such a neutralino 
can give an appropriate contribution to the dark matter density consistent with the recent 
observations if it has mass 35 — 50 GeV [23|. In this case the lightest Higgs boson decays 
predominantly into Inert neutralino whereas usual Higgs branching ratios are less than a 
few percent. 

The layout of this paper is as follows. In Section 2 we briefly review the EgSSM. In 
Sections 3 and 4 the spectrum and couplings of the Inert neutralinos, charginos and Higgs 
bosons are specified. The novel decays of the lightest CP-even Higgs state and dark matter 
constraints are discussed in Section 5. In section 6 we discuss the benchmark points and 
the experimental constraints and predictions. Section 7 summarizes and concludes the 
paper. 



2. Exceptional SUSY model 

The EeSSM is based on the SU{3)c x SU{2)w x f/(l)y x U{1)n gauge group which is a 
subgroup of Eq. The additional low energy U{1)n, which is not present either in the SM 
or in the MSSM, is a linear superposition of U{1)^ and U{1)^, namely 

U{1)n = \u{1)^ + ^U{1)^, (1) 

where the U{1)^ and U{1)^ symmetries are defined by: 

Ee ^O(IO) X U{1)^ , 50(10) ^ SU{^) x f/(l)^ . 

Thus the EgSSM can originate from an i^g GUT gauge group which is broken at the GUT 
scale Mx- The extra U{1)n gauge symmetry is defined such that right-handed neutrinos 
carry zero charges. 

In Eq theories the anomalies cancel automatically; all models that are based on the 
Eq subgroups and contain complete representations of Eq should be anomaly-free. Con- 
sequently, in order to make a supersymmetric model with an extra U{1)n anomaly-free, 
one is forced to augment the minimal particle spectrum by a number of exotics which, 
together with ordinary quarks and leptons, form complete fundamental 27 representations 
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27i ^ 10, ^ + 5*, ^ + 5*, = + 5, ^ 



of £^6- Thus the particle content of the EgSSM involves at least three fundamental repre- 
sentations of Eq at low energies. These multiplets decompose under the SU{5) x U{1)n 
subgroup of Eq as follows [27] : 

(2) 

, ] +(l,Ol. . 

The first and second quantities in brackets are the SU{5) representation and extra U{1)n 
charge respectively, while z is a family index that runs from 1 to 3. An ordinary SM family, 
which contains the doublets of left-handed quarks Qi and leptons Li, right-handed up- 
and down-quarks {u^ and d^) as well as right-handed charged leptons, is assigned to 



10, —= + 5*, —= . Right-handed neutrinos should be associated with the 



last term in Eq. (|2]), (1, 0)^. Because they do not carry any charges right-handed neutrinos 

are expected to be superheavy allowing them to be used for both the see-saw mechanism 
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and leptogenesis. The next-to-last term, ( 1, — = ) , represents SM-singlet fields Si. 



which carry non-zero U{l)]\f charges and therefore survive down to the EW scale. The pair 
of S'f/(2)iy-doublets {Hf and H^) that are contained in [5*, = J and [5 



/40/i V V40, 

have the quantum numbers of Higgs doublets. They form either Higgs or Inert Higgs 
SU{2)iv multiplets §. Other components of these SU{5) multiplets form colour triplets 
of exotic quarks Di and Di with electric charges —1/3 and +1/3, respectively. These 
exotic quark states carry a B — L charge ±2/3, twice that of ordinary ones. Therefore in 
phenomenologically viable Eq inspired models they can be either diquarks or leptoquarks. 

In addition to the complete 27, multiplets the low energy matter content of the EgSSM 
can be supplemented by an SU{2)w doublet I/4 and anti-doublet L4 from the extra 27' 
and 27' to preserve gauge coupling unification. These components of the Eq fundamental 
representation originate from (5*, , | of 27' and | 5, = | of 27' by construction. 

V V40/ V V4o; 

Thus, in addition to a Z' corresponding to the U{1)n symmetry, the EgSSM involves 
extra matter beyond the MSSM that forms three 5 + 5* representations of SU{5) plus 
three SU{5) singlets with U{1)n charges. The analysis performed in [2B] shows that the 
unification of gauge couplings in the EgSSM can be achieved for any phenomenologically 
acceptable value of 0:3 (M^) consistent with the measured low energy central value, unlike 
in the MSSM which, ignoring the effects of high energy threshold corrections, requires val- 
ues of 03 (M^) which are significantly above the experimentally measured central value. 
The presence of a Z' boson and of exotic quarks predicted by the EgSSM provides spec- 
tacular new physics signals at the LHC which were discussed in [23]-[2l], [29]. Recently 



^We use the terminology "Inert Higgs" to denote Higgs-like doublets that do not develop VEVs. 
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the particle spectrum and collider signatures associated with it were studied within the 
constrained version of the EgSSM |30] . 

In general, the Eq symmetry does not forbid lepton and baryon number violating 
operators that result in rapid proton decay. Moreover, exotic particles in Eq inspired 
SUSY models give rise to new Yukawa interactions that induce unacceptably large non- 
diagonal flavour transitions. To suppress these effects in the EgSSM an approximate Z2 
symmetry is imposed. Under this symmetry all superfields except one pair of Hf and if" 
(say = and = H^) and one SM-type singlet field {S = S3) are odd. The Z2 
symmetry reduces the structure of the Yukawa interactions to 

+ ^M,,N^N;; + fi\L,i,) + hf^{HMe', + h^,{HuL^^ (3) 

where a, (3 = 1,2 and i,j = 1,2,3. The SU{2)w doublets Hu and Hd and SM-type 
singlet field 5*, that are even under the Z2 symmetry, play the role of Higgs fields. At 
the physical vacuum the Higgs fields develop VEVs 

generating the masses of the quarks and leptons. Instead of Vi and f 2 it is more convenient 





to use tan/3 = V2/V1 and v = \fv\^rv\ = 246 GeV. The VEV of the SM-type singlet 
field, s, breaks the extra U{1)n symmetry thereby providing an effective /i term as well 
as the necessary exotic fermion masses and also inducing that of the Z' boson. Therefore 
the singlet field 5* must acquire a large VEV in order to avoid conflict with direct particle 
searches at present and past accelerators. This also requires the Yukawa couplings Aj 
and Ki to be reasonably large. If Aj or Ki are large enough at the GUT scale they affect 
the evolution of the soft scalar mass m| of the singlet field S rather strongly resulting 
in a negative value of m| at low energies which triggers the breakdown of the U{1)n 
symmetry. 

Note that the surviving components from the 27' and 27' manifest themselves in the 
Yukawa interactions ([3]) as fields with lepton number L = ±1. The corresponding mass 
term ^'L^^L^^ in the superpotential ([3]) is not involved in the process of electroweak sym- 
metry breaking (EWSB). Moreover this term is not suppressed by the Eq symmetry. 
Therefore the parameter /i' remains arbitrary. Gauge coupling unification requires ^' to 
be below about 100 TeV |28] . Thus we assume that the scalar and fermion components 
of the superfields I/4 and L4 are very heavy so that they decouple from the rest of the 
particle spectrum. 
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Although Z2 ehminates any problems related with baryon number violation and non- 
diagonal flavour transitions it also forbids all Yukawa interactions that would allow the 
exotic quarks to decay. Since models with stable charged exotic particles are ruled out 
by various experiments [21] the symmetry must be broken. At the same time, the 
breakdown of should not give rise to operators that would lead to rapid proton decay. 
There are two ways to overcome this problem: the Lagrangian must be invariant with re- 
spect to either a Z^ symmetry, under which all superflelds except leptons are even (Model 
I), or a Z^ discrete symmetry, which implies that exotic quark and lepton superflelds are 
odd whereas the others remain even (Model II). If the Lagrangian is invariant under the 
Z^ symmetry, then the terms in the superpotential which permit exotic quarks to decay 
and are allowed by the symmetry can be written in the form 



that implies that exotic quarks are diquarks. If Z2 is imposed then the following couplings 
are allowed: 



In this case baryon number conservation requires the exotic quarks to be leptoquarks. 

3. Inert charginos and neutralinos 

From here on we assume that Z2 symmetry violating couplings are small and can be 
neglected in our analysis. This assumption can be justifled if we take into account that the 
Z2 symmetry violating operators may give an appreciable contribution to the amplitude 
of — oscillations and give rise to new muon decay channels like /i — )■ e~e"'"e~. In 
order to suppress processes with non-diagonal flavour transitions the Yukawa couplings 
of the exotic particles to the quarks and leptons of the flrst two generations should be 
smaller than 10^^ — 10^^. Such small Z2 symmetry violating couplings can be ignored 
in the flrst approximation. 

In this approximation and given the previous assumption that only Hu-, Hd and S ac- 
quire non-zero VEVs the charged components of the Inert Higgsinos (-^2^, H"^^ , -^2" 5 
and ordinary chargino states do not mix. The neutral components of the Inert Higgsinos 
{Hf^^ Hf^ , H^^, -^2°) Inert singlinos {Si, S2) also do not mix with the ordinary 
neutralino states. Moreover if Z2 symmetry was exact then both the lightest state in the 
ordinary neutralino sector and the lightest Inert neutralino would be absolutely stable. 
Therefore, although Z2 symmetry violating couplings are expected to be rather small, 
we shall assume that they are large enough to allow either the lightest neutralino state 




A(qA) + 4.Ac^>l 



(5) 



W2 = 9f,,etD,ul + gg,iQ,L,)D,. 



(6) 
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or the lightest Inert neutrahno to decay within a reasonable time, the lighter of the two 
being the stable LSP. 

In the field basis (i^f'^, i?2°! '^2, Hf^, H'f, Si) the mass matrix of the Inert neutralino 
sector takes a form 

MiN=\ , (7) 

V ^12 ^11 / 

where A^p are 3x3 sub-matrices given by [26] : 

/ Xaps fi3aV sin (3 \ 

A/3aS fi3aV COS (3 

\ fafivsmf3 fapv cos 15 / 

so that Aap = A^^. In the basis of Inert chargino interaction states {H2~^, , , Hf 
the corresponding mass matrix can be written as 



- "7! 



a 



1 



Mic=\ , C^p = s . (9) 



C 



V2 



where Cap are 2x2 sub-matrices. From Eqs. ([7j)-(j9]) one can see that in the exact Z2 
symmetry limit the spectrum of the Inert neutralinos and charginos in the EgSSM can be 
parametrised in terms of 

Aa/3, /a/3, /a/3, tan/3, S. (10) 

In other words the masses and couplings of the Inert neutralinos are determined by 12 
Yukawa couplings, which can be complex, tan /3 and s. Four of the Yukawa couplings 
mentioned above, i.e. A^^, as well as the VEV of the SM singlet field s set the masses 
and couplings of the Inert chargino states. Six off-diagonal Yukawa couplings define the 
mixing between the two families of the Inert Higgsinos and singlinos. 

In the following analysis we shall choose the VEV of the SM singlet field to be 
large enough (s > 2400 GeV) so that the experimental constraints on Z' boson mass 
{Mz' ^ 865 GeV) and Z — Z' mixing are satisfied. In order to avoid the LEP lower limit 
on the masses of Inert charginos the Yukawa couplings Xap are chosen so that all Inert 
chargino states are heavier than 100 GeV. In addition, we also require the validity of 
perturbation theory up to the GUT scale and that constrains the allowed range of all 
Yukawa couplings. 

The theoretical and experimental restrictions specified above set very strong limits on 
the masses and couplings of the lightest Inert neutralinos. In particular, our numerical 
analysis indicates that the lightest and second hghtest Inert neutralinos are always light. 
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They typically have masses below 60 — 65 GeV. These neutralinos are predominantly Inert 
singlinos. From our numerical analysis it follows that the lightest and second lightest Inert 
neutralinos might have rather small couplings to the Z-boson so that any possible signal 
which these neutralinos could give rise to at LEP would be extremely suppressed. As a 
consequence such Inert neutralinos would remain undetected. At the same time four other 
Inert neutralinos, which are approximately linear superpositions of neutral components 
of Inert Higgsinos, are normally heavier than 100 GeV. 



3.1 The diagonal Inert Yukawa approximation 



In order to clarify the results of our numerical analysis, it is useful to consider a few simple 
cases that give some analytical understanding of our calculations. In the simplest case 
when all off-diagonal Yukawa couplings vanish, considered in 



A. 



Ao 5a/3, 



/a/3 — fa 5al3i faP — fa 



the mass matrix of Inert neutralinos reduces to the block diagonal form while the masses 
of the Inert charginos are given by 

m ± = . (11) 

In the limit where fa = fa one can easily prove using the method proposed in [32] 
that there are theoretical upper bounds on the masses of the lightest and second lightest 
Inert neutralino states. The corresponding theoretical restrictions are 



I Xa I 1^ 



/^a=2 



r2 2 

|m^±|2 + :^(l + sin2 2/3) ) - fy sm^2P 



(12) 



The value of /i^ decreases with increasing \fn^±\ and tan/3, hence reaching its maximum 
value of -^v for m ± and tan/3 1. At large values of {m ±\ and tan /3, Eq. 0121) 
simplifies resulting in 



I Xa I ~ 



fy sin' 2/3 



m^±P + ^(l + sin2 2/3) 



(13) 



Eqs. fll2l) - fll3l) demonstrate that the upper bound on the mass of the lightest Inert neu- 
tralino also depends on the values of Yukawa couplings and fa- For relatively small 
values of tan/3, the theoretical restrictions on fa and fa, due to the requirement that the 
perturbation theory is valid up to the GUT scale, become weaker with increasing tan /3. 
However, at large values of tan f3 the upper bounds on |m^o | become rather small according 
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to Eqs. f|T2l) - f|T3l) . When tan (3 tends to unity, /i^ also decreases because the constraints on 
fa and fa become more and more stringent. The theoretical restrictions on \Tn^o\ achieve 
their maximal value around tan /3 ^ 1.5. For this value of tan /3 the requirement of the 
validity of perturbation theory up to the GUT scale implies that /i = /i = /2 = /2 are 
less than 0.6. As a consequence the lightest Inert neutralinos are lighter than 60 — 65 GeV 
for |m +1 > 100 GeV. 

The Inert neutralino mass matrix ([7])- ([8]) can be diagonalized using the neutralino 
mixing matrix defined by 



no sum on i. 



(14) 



In the limit where off-diagonal Yukawa couplings vanish and A^s ^ faV, faV the eigen- 
values of the Inert neutralino mass matrix can be easily calculated (see [26]). The masses 
of the four heaviest Inert neutralinos are set by the masses of Inert chargino states 

fafav'^ sin 2 (3 



A.3,4,5,6 



4:171 + 

Xa 



(15) 



The masses of the two lightest Inert neutralinos are determined by the values of the 
Yukawa couplings and 

fafav'^sm2(3 



A.OL 



2m ± 

Xa 



(16) 



These are naturally small and hence good candidates for being the LSP and NLSP since 



m ± 

Xa 



s from Eq. fill I) and hence m^o ~ v'^/s as observed in 



Again one can see that the masses of the lightest Inert neutralino states decrease with 
increasing tan /3 and chargino masses. In this approximation the lightest Inert neutralinos 
are made up of the following superposition of interaction states 



Aa 



NlHf + NlHf + NIS2 + KHf + NlHf + NlSi , 



(17) 



where 



AT 1 = Nf = Nf = 0. 



~ 1 



^ 2 \\is 

f2V COS /3 



2 Wos 



fiv cos /3 



fl cos^ /3 + fl sin^ /3 
f2V sin /3 



^1 



A2S 



/|cos2/3 + /|sin2/3 



/if sin /3 
Ais 



iV| = iV| = iV| = 0, 



From Eq. (fT8|) it becomes clear that the lightest and second lightest Inert neutralinos are 
mostly Inert singlinos. 
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Using the above lightest and second hghtest Inert neutrahno compositions it is straight- 
forward to derive the couphngs of these states to the Z-boson. In general the part of the 
Lagrangian that describes the interactions of Z with x'j* and can be presented in the 
following form: 



(19) 



In the case where off-diagonal Yukawa couplings go to zero while XaS ^ faV, faV the 
relative couplings of the lightest and second lightest Inert neutralino states to the Z-boson 
are given by 



Rzaf} = Rzaa ^afi , Rzaa = ( fa ^^S^ - fl sin^ /3 ) . (20) 



Eq. ( 120|) demonstrates that the couplings of x? and X2 the Z-boson can be very strongly 
suppressed or even tend to zero. This happens when \fa \ cos/3 = sin/3, which is when 
Xq contains a completely symmetric combination of H'^ and H'^. Eq. (120 p also indicates 
that the couplings of x? and X2 Z are always small when Inert charginos are rather 
heavy or fa and fa are small (i.e. m-^o — )■ 0). 

3.2 A27 and pseudo-Dirac lightest neutralino states 

In order to provide an explanation of the origin of the aproximate Z2 symmetry that sin- 
gles out the third family of Higgs doublets and singlets, and to account for tri-bimaximal 
mixing and other features of the quark and lepton spectrum, a A27 family symmetry has 
been applied to the EgSSM [33] 0. The addition of a A27 family symmetry implies an Inert 
neutralino mass matrix with An ^ A22 ~ 0, where Aap are defined in Eq. ([7]), leading to 
approximately degenerate lightest neutralino states with a pseudo-Dirac structure. 

When all flavour diagonal Yukawa couplings A^a, faa and exactly vanish, i.e. 
All = ^22 = 0, all Inert Higgsinos and singlinos form Dirac states. In this limit the 
Lagrangian of the EgSSM is invariant under a U{1) global symmetry. The fermion com- 
ponents of the Inert Higgs superfields transform under this symmetry as follows: 

Si e'^Si , e'^^H^ , Hf e^'^Hf , 



"^The corresponding mass terms come from the product (3 x 3 x 3/)(3 x 3^), where 3, 3/ and 3^ are 
triplet representations of A27. The 21 i multiplets that contain quarks and leptons form 3 representations 
of the A27 group. 3/ and 3^ contain flavon fields that break A27. In the considered model the non-zero 
mass of the lightest Inert neutralino state is induced by the symmetric invariant that appears in the 
(3x3x3) decomposition of the A27 triplet representation (i.e. 123 + 231 + 312 + 213 + 321 + 132) [31]. 
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In the above limiting case the hghtest Inert neutralino is a Dirac state formed predom- 
inantly by 5*1 and 5*2- In this case the LSP and its antiparticle have opposite charges 
with respect to the extra global U{1) and this might lead to so-called asymmetric dark 
matter (ADM) In the framework of the ADM scenario there can be an asym- 

metry between the density of dark matter particles and their antiparticles in the early 
universe similar to that for ordinary baryons. This may have a considerable effect on the 
relic density calculations [36]. In particular, if an asymmetry exists between the number 
density of dark matter particles and their antiparticles in the early universe, then one 
can get an appreciable dark matter density even if the dark matter particle-antiparticle 
annihilation cross section is very large like in the case of baryons. Moreover if most of the 
dark matter antiparticles are eliminated by annihilation with their particles then such an 
ADM scenario does not have the usual indirect signatures associated with the presence 
of dark matter (e.g. there is no high energy neutrino signal from annihilations in the Sun 
etc.). At the same time, a relatively high concentration of dark matter particles can build 
up in the Sun altering heat transport in the solar interior and affecting the low energy 
neutrino fluxes [57] . 

In practice the A27 scenario tells us that we are somewhat away from the above limit- 
ing case, with a broken global U{1) symmetry leading to almost degenerate pseudo-Dirac 
lightest neutralinos, where the relic density of the LSP can be calculated by standard 
methods. It will turn out that the LSP cannot be too light (must be of order M^/2) 
in order not to have too high a cosmological relic density. At the same time we will see 
that the two lightest neutralinos cannot be too heavy in order for perturbation theory 
to be valid up to the GUT scale. In practice this means that in realistic scenarios the 
two lightest Inert neutralino states are rather close in mass. The A27 scenario provides a 
natural explanation of this successful neutralino mass pattern. It is worth noting that the 
results from the previous section can be reinterpreted in terms of this scenario. Specif- 
ically in the case where An = A22 = and A21 = Ayi a block diagonalisation of the 
Inert neutralino mass matrix ([7]) yields ^22 ~ and A^i — )■ = A21 (with 

A21 = A\2 — ^21 = ^'12 — 0)- This only corresponds to a redefinition of the generations 1 
and 2 and does not mix fields of different hypercharge. This provides the following dictio- 
nary between these two scenarios: A'^ = -A22 = A21; /n = -/22 = /21; /n = -/22 = h\- 
Rewriting the Inert neutralino mass matrix in this block diagonal form also makes it clear 
that the Rzn coupling vanishes in this limit, as it did in the subsection (13. ip . 

3.3 Scenario with one light family of Inert Higgsinos 

Another limit that it is worth considering corresponds to the case where one pair of Inert 
Higgs doublets decouples from the rest of the spectrum. This occurs when either the 
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corresponding states are extremely heavy (> 1 TeV) or they have rather small couplings 
to other Inert Higgs fields. When ifg" ifg ° decouple, the Inert neutralino mass matrix 
([7]) reduces to a 4 x 4 matrix. If Ans ^ faiv, faiv, the Inert Higgs states associated with 
Hf^ and can be integrated out. Then the resulting 2x2 mass matrix can be written 
as follows 

2/11/11 /11/21 + /11/21 

/11/21 + /11/21 2/21/21 
The masses of the lightest and second lightest Inert neutralinos, which are predominantly 
superpositions of the Inert singlinos 5*1 and S2, are given by 



IS 



v'^ sin 2/3 
4m, ± 

Xi 



(22) 



^x?,xg 



sin 2/3 

4:771 + 

Al 



/ii/ii + /21/21 ± V(/ii + /li)(/ii + /li) 



(23) 



From Eq. (123|) it is easy to see that the substantial masses of the lightest and second 
lightest Inert neutralinos can be induced even if only one family of the Inert Higgsinos 
couples to Si and 82- 

Using Eq. (IT^ one can also calculate the couplings of Xi and X2 ^^e Z-boson 
,2 r 

(/ii cos 9 + /21 sin cos^ /3 - (/n cos 9 + /21 sin 6^ sin^ /3 



Rzii 

RZ22 
Rzi2 



2m2 



2m2 



xt 



(/21 cos^ — /ii sin6')^cos^/3 — (/21COS6' — /n sin 5)^ sin^ /3 



R 



Z2l 



2m2± 

Xi 



(24) 



2(/li - /n) sin2^ + /n/21 cos2^ ) cos^ (3 
- (lifii - III) sin2^^ + /n/21 cos2e ) sin^/3 



where tan 26* = (/11/21 + /2i/ii)/(/2i/2i - /ii/ii)- Again from Eqs. ([24]) it follows that 
-Rzii, -Rz22 and Rzi2 are typically small since m^± ~ s from Eq. (ITTi) and hence they 
are proportional to v'^/s^. However this assumes the lightest Inert chargino is rather 
heavy. Alternatively the couplings may be small due to a cancellation between different 
contributions in Eqs. flM|l . and/or the f-couplings being small (i.e. ^x1,X2 ~^ 

The simple hierarchical structure of the spectrum of the Inert neutralinos considered 
above allows us to highlight an interesting scenario which does not normally appear in the 
simplest SUSY extensions of the SM such as the MSSM and NMSSM. When /n = /21 = 
the diagonal entries of the mass matrix (122|) vanish leading to the formation of a Dirac 
lightest Inert neutralino state. In this case the Lagrangian of the model is invariant under 
extra U{1) global symmetry transformations [f| Si — )■ e^'^S'i , Hi — )■ e~*"-ff" , Hf — )■ e^°'Hf , 
^2 — )■ e~*"S'2 . In fact if the EgSSM possess such an exact U{1) global symmetry, then the 
spectrum of the Inert neutralinos contains a set of Dirac states only. 



^Similar results can be obtained for /n = /21 = 
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4. Higgs masses and couplings 



The presence of light Inert neutrahnos in the particle spectrum of the EgSSM makes 
possible the decays of the Higgs bosons into these exotic final states. In this and the next 
section we argue that such decays may result in the modification of the SM-hke Higgs 
signal at current and future colliders. Since our main concern in this paper is the decays of 
the SM-like lightest Higgs boson, we shall ignore the effects of the Inert Higgs scalars and 
pseudoscalars which do not mix appreciably with the scalar sector responsible for EWSB. 
We also assume that all the Inert bosons are heavier than the SM-hke Higgs boson. 

The sector responsible for the EWSB in the EgSSM includes two Higgs doublets Hy, 
and Hd as well as the SM singlet field S. The Higgs effective potential can be written in 
the following form: 

V = Vp + VD + Vsoft + ^V , 



'■uj 

,2 . .9 „/2 



(25) 



Vsoft = ml\S\'^ + ml\Hd\'^ + ml\Hy,\'^ + 



XAxS{HuHd) + h.c. 



where §2, g' = a/3/5(7i and g[ are the low energy SU{2)\y, U{1)y and f/(l)jv gauge 
couplings while Qi, Q2 and Qs are the effective U{1)n charges of Hd, and S. The 
term AV represents the contribution from loop corrections to the Higgs effective potential. 
Here = {H^ H'), = (H^, H'J and (//,//J = H^H- - H'^H^ 

Initially the EWSB sector involves ten degrees of freedom. However four of them are 
massless Goldstone modes which are swallowed by the W^, Z and Z' gauge bosons that 
gain non-zero masses. In the hmit where s ^ v the masses of the W^, Z and Z' gauge 
bosons are given by 



where g = \/gf+g^. When CP-invariance is preserved the other degrees of freedom form 
two charged, one CP-odd and three CP-even Higgs states. The masses of the charged 
and CP-odd Higgs bosons are 



- - + + A„ < . --^s + A. , (26) 

where A± and are the loop corrections. 
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The CP-even Higgs sector involves ReH^, ReH^ and ReS. In the field space basis 
[h, H, N), rotated by an angle /3 with respect to the initial one, 



ReH^ = {hcosl3 - Hsml3 + vi)/V2, 
ReHl = {hsm/3 + Hcos(3 + V2)/V2, 
ReS = {s + N)/V2, 
the mass matrix of the CP-even Higgs sector takes the form [38] : 
/ d^V 1 d^V d^V \ 



where 



Ml 



1 d^V 



V dvdf3 dvds 
1 d^V 1 d^V 



V dvd/S f 2 (9^/3 V dsd/S 
d^V 1 d^V d^V 



( 



Mt, Mt^ Mt, 



2 \ 



Ml Ml 
Ml Ml 



Ml 



M? 



33 / 



\ dvds v dsdp d'^s J 



2 



(27) 



(28) 



'—v"^ sin^ 2/3 + cos^ 2/3 + g^v'^{Qi cos^ /3 + Q2 sin^ pf + An , 



Ml 



4 

4 8/ ^ 2 



gi)x 



X (Qi cos^ I3 + Q2 sin^ /3) sin 2/3 + A12 , 
^h^^+U"^J' 2/3+ -(Q, 



gov sin' 2/3 + A22 



Ml 



Ml 



Ml 



Ml 



Ml 



COS 2/3 + ^(^52 



gi)Q5t;ssin2/3 + A23 



■V sin 2/3 + A^f s + g[^{Qi cos^ /3 + sin^ l3)Qsvs + A13 , (29) 



^t;2sin2/3 + M|, + A33 
2y/2s 



In Eq. ( 129|) the A^- represent the contributions from loop corrections which in the leading 
one-loop approximation are rather similar to the ones calculated in the NMSSM@. bmce 



the minimal eigenvalue of the mass matrix f l28p -( 12^ is always less than its smallest 
diagonal element, at least one Higgs scalar in the CP-even sector (approximately h) 
remains always light, i.e. m^^ < Ml. In the leading two-loop approximation the mass of 
the lightest Higgs boson in the EgSSM does not exceed 150 — 155 GeV. When the SUSY 
breaking scale Ms and the VEV s of the singlet field are considerably larger than the EW 
scale, the mass matrix fl2S]) - fl2^ has a hierarchical structure and can be diagonalised using 



®Note that the exphcit expressions for Ay, A± and presented in the first paper in [38j contain a 
typo. In the corresponding formulae /i is neither a parameter of the MSSM Lagrangian nor an effective 
/Lt-term in the NMSSM. It has to be associated with the renormalisation scale. 
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the perturbation theory [38]-|in]- In this case the masses of the heaviest Higgs bosons are 
closely approximated by the diagonal entries Mfg and MI3 [23]. As a result the mass of 
one CP-even Higgs boson (approximately given by H) is governed by while the mass 
of another one (predominantly the N singlet field) is set by Mz'- When A > (7^, vacuum 
stability requires rriA to be considerably larger than Mz> and the EW scale so that the 
qualitative pattern of the Higgs spectrum is rather similar to the one which arises in the 
PQ symmetric NMSSM [1Q]-[1T]. In the considered limit the heaviest CP-even, CP-odd 
and charged states are almost degenerate around and lie beyond the TeV range [25] . 

If all other Higgs states are much heavier than the lightest CP-even Higgs boson then 
the lightest Higgs state (approximately given by h) manifests itself in the interactions 
with gauge bosons and fermions as a SM-like Higgs boson. Since within the EgSSM the 
mass of this state is predicted to be relatively low its production cross section at the LHC 
should be large enough so that it can be observed in the near future. In this context it is 
particularly interesting and important to analyse the decay modes of the lightest CP-even 
Higgs state. Furthermore we concentrate on the decays of the SM-like Higgs boson into 
the lightest and second lightest Inert neutralinos. 

The couplings of the Higgs states to the Inert neutralinos originate from the inter- 
actions of Hu, Hd and 5* with the Inert Higgs superfields in the superpotential. Using 
Eqs. fl27j) one can express ReH^, ReH^ and Re S in terms of the components of the 
CP-even Higgs basis h, H and A^. At the same time the components of the CP-even 
Higgs basis are related to the physical CP-even Higgs eigenstates by virtue of a unitary 
transformation: 



Combining all these expressions together one obtains an effective Lagrangian that de- 
scribes the interactions of the Inert neutralinos with the CP-even Higgs eigenstates 



/ h \ 



(hA 



H 



(30) 





(31) 



= hiNfN^ + h2NfN] + f2iNfN^ + f22NfN] , 
4 = /nATf AT/ + fuNfN] + hiNfNj + h^NfN] , 
A,,- = XnNfN^ + + \2iN}N^ + \22N}N] , 
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where i,j = 1,2, ...6 and m = 1,2,3. In Eq. fl3T]) = {—i'j^Y'x'i is the set of Inert 
neutralino eigenstates with positive eigenvalues, while 9i equals (1) if the eigenvalue 
corresponding to Xi is positive (negative). As before, the Inert neutralinos are labeled 
according to increasing absolute value of mass, with ijj^ being the lightest Inert neutralino 
and ipe the heaviest. 

The expressions for the couplings of the Higgs scalars to the Inert neutralinos f l3ip 
become much more simple in the case of the hierarchical structure of the Higgs spectrum. 
In this case Uij is almost an identity matrix. As a consequence, the couplings of the SM-like 
Higgs boson to the lightest and second lightest Inert neutralino states are approximately 
given by 



(32) 



where 7, cr = 1,2, labeling the two light, mostly Inert singlino states. In the limit when 
off-diagonal Yukawa couplings that determine the interactions of the inert Higgs fields 
with Hu, Hd and S vanish, as defined in subsection fl3.ip . and Inert neutralino mass 
matrix has a hierarchical structure (i.e. A^s ^ /af , /of), one can use the expressions 
f lTSj) for in order to derive the approximate analytical formulae for X^^. Substituting 
Eqs. f lTS]) into f l5^ one obtains 



X^^l ^ ^ 5,. , (33) 

These simple analytical expressions for the couplings of the SM-like Higgs boson to the 
lightest and second lightest Inert neutralinos are not as surprising as they may first ap- 
pear. When the Higgs spectrum is hierarchical, the VEV of the lightest CP-even state 
is responsible for all light fermion masses in the EgSSM. As a result we expect that their 
couplings to SM-like Higgs can be written as usual as being proportional to the mass 
divided by the VEV. We see that this is exactly what is found in the limit of \ being 
small. 



5. Novel Higgs decays and Dark Matter 
5.1 Higgs decay widths 

The interaction Lagrangian ( 13T|) gives rise to decays of the lightest Higgs boson into Inert 
neutralino pairs with partial widths given by 

2 

' (34) 



I 19 1^ |2n 2 I™ |2|™ 12 
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where = - (1) for a = /3 (a 7^ 

The partial widths associated with the exotic decays of the SM-hke Higgs boson fl34l) 
have to be compared with the Higgs decay rates into the SM particles. When the SM-like 
Higgs state is relatively light (m/j^ < 140 GeV) it decays predominantly into 6-quark and 
r-lepton pairs. The partial decay width of the lightest CP-even Higgs boson into fermion 
pairs is given by (for recent review see [32]) 



nhi -> //) = 



m 



9l 



327r V M, 



/ 



1 



4m^ 



3/2 



(35) 



m/ \ ""hi 

Eq. ( 135|) can be used for the calculation of the lightest Higgs decay rate into r-lepton pairs. 
In this case the coupling of the lightest CP-even Higgs state to the r-lepton normalized 
to the corresponding SM coupling, i.e. ghirr, is given by 

1 



9hn 



cos/3 



Ulh, COS sin/3 



(36) 



For a final state that involves 6-quarks one has to include the QCD corrections. In 
particular, the fermion mass in Eq. ( 135|) should be associated with the running 6-quark 
mass mb(yu). The bulk of the QCD corrections are absorbed by using the running 6-quark 
mass defined at the appropriate renormalisation scale, i.e. at the scale of the lightest 
Higgs boson mass {fi = rrihi) ^^e considered case. In addition to the corrections which 
are associated with the running 6-quark mass there are other QCD corrections to the 
Higgs coupling to the 6-quark that should be taken into account [43j. As a consequence, 
the partial decay width of the lightest CP-even Higgs boson into 6-quark pairs can be 
calculated using Eq. (135|) if one sets Nc = 3 and replaces 



rrif mb[mh^^ 
1 



cos^ /3 



f^Lcos/3-f/]^^^ sin/3 



H 



A 



bb 



5.67— + (35.94 - 1.36iV.)% 



(37) 



Ah^^[ 1.57- -loj 



where a. 



^h, 1,2 "^b 

^ + -log^^ 



m 



m 



hi 



Here we neglect radiative corrections that originate from loop 



diagrams that contain SUSY and exotic particles 0. 

From Eqs. f l33p - (|5S]) one can see that in the EgSSM the branching ratios of the SM-like 
Higgs state into the lightest and second lightest Inert neutralinos depend rather strongly 
on the masses of these exotic particles. When the lightest Inert neutralino states are 



^Radiative corrections that are induced by SUSY particles can be very important particularly in the 
case of the bottom quark at high values of tan/3 (for a review, see IH]). 
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relatively heavy, i.e. 'mxi,X2 ^ ^bi^hj, the lightest Higgs boson decays predominantly 
into XaXi3 while the branching ratios for decays into SM particles are suppressed. On the 
other hand if the lightest Inert neutralinos have masses which are considerably smaller 
than the masses of the 6-quark and r-lepton then the branching ratios of the exotic 
decays of the SM-like Higgs state are small. In the EgSSM the lightest and second 
lightest Inert neutralinos are expected to be heavier than a few MeV so that they would 
not contribute to the expansion rate prior to nucleosynthesis and thus not modify Big 
Bang nucleosynthesis (BBN). 

5.2 Dark matter 

More stringent constraints on the masses of the lightest Inert neutralino can be obtained 
if we require that this exotic state accounts for all or some of the observed dark matter 
relic density which is measured to be f2cDM^^ = 0.1099 ± 0.0062 [45]. If a theory predicts 
a greater relic density of dark matter than this then it is ruled out, assuming standard 
pre-BBN cosmology. A theory that predicts less dark matter cannot be ruled out in the 
same way but then there would have to be other contributions to the dark matter relic 
density. 

In the limit where all non-SM fields other than the two lightest Inert neutralinos are 
heavy (> TeV) the lightest Inert neutralino state in the EgSSM results in too large a 
density of dark matter. As we noted in Section 3, Xi is usually composed of Inert singlino 
and has a mass (Eq. (16)) which is inversely proportional to the charged Higgsino mass. 
Thus in this limit it is typically very light | m^o \ -C Mz ■ As a result the couplings of the 
lightest Inert neutralino to gauge bosons, the SM-like Higgs state, quarks and leptons are 
quite small leading to a relatively small annihilation cross section for XiXi ~^ SM particles. 
Since the dark matter number density is inversely proportional to the annihilation cross 
section at the freeze-out temperature (see, for example [l6j) the lightest Inert neutralino 
state gives rise to a relic density that is typically much larger than its measured value. 
Thus in the limit considered the bulk of the EgSSM parameter space that leads to small 
masses of Xi is ruled out. 

The situation changes dramatically when the mass of the lightest Inert neutralino 
increases. In this case the Higgsino components of Xi become larger and as a consequence 
the couplings of Xi to the Z-boson grow [26j. A reasonable density of dark matter can 
be obtained for \itl^o \ ~ Mz/2 when the lightest Inert neutralino states annihilate mainly 
through an s-channel Z-boson, via its Inert Higgsino doublet components which couple 
to the Z-boson. It is worth noting that if Xi was pure Inert Higgsino then the s-channel 
Z-boson annihilation would proceed with the full gauge coupling strength leaving the 
relic density too low to account for the observed dark matter. In the EgSSM the LSP is 
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Figure 1: Contour plot of {X'^^)'^ and relic density ^yh? regions in the {f , tan /3)-plane with 
s = 2400 GeV, = = A„„ = 0, fu = f, fu = fu/a, /21 = 1.02 • /12, hi = 0.98 • fu, 
a = 0.75 + 0.25 tan /3 and A12 = A21 = 0.06 (^^iT-^i = 101.8 GeV^. The red region is where the 
prediction for VL^K^ is consistent with the measured one a range o/J^cdm/i^ = 0.1099 it 0.0062. 
The dark green region corresponds to D < 3 (D is defined in subsection 6.2) while the pale green 
region represents the part of the parameter space in which D varies from 3 to 4. The grey area 
indicates that D > A. The blue region corresponds to m^o > Mzjl, while the dark blue region to 
the right is ruled out by the requirement that perturbation theory remains valid up to the GUT 
scale. 

mostly Inert singlino so that its coupling to the Z-boson is typically suppressed, since it 
only couples through its Inert Higgsino admixture leading to an increased relic density. 
In practice, the appropriate value of ^cdm^^ can be achieved even if the coupling of x? 
to the Z-boson is relatively small. This happens when x5 annihilation proceeds through 
the Z-boson resonance, i.e. 2|m^o| ^ Mz [SHI HZ]. Thus scenarios which result in a 
reasonable dark matter density correspond to lightest Inert neutralino masses that are 
much larger than mh{mh^), and hence the SM-like Higgs has very small branching ratios 
into SM particles. 



6. Benchmarks, constraints and predictions 

In order to illustrate the features of the EgSSM mentioned in the previous section, we 
shall specify a set of benchmark points (see Tables 1-2). For each benchmark scenario 
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Figure 2: Contour plot of {X'^^)'^ and relic density ^yh? regions in the {f,tan /3)-plane with 
s = 2400 GeV, = = A„„ = 0, fu = f, fi2 = fu/a, hi = 1-02 • /12, hi = 0.98 • /12, 
a = 0.5 + 0.5 tan /3 and A12 = A21 = 0.06 fw-^± = 101.8 GeVj. The red region is where the 
prediction for i^^h'^ is consistent with the measured one a range o/ricDM^i^ = 0.1099 it 0.0062. 
The dark green region corresponds to D < 3 while the pale green region represents the part of the 
parameter space in which D varies from 3 to 4. The blue region corresponds to m^o > Mz/2, 
while the dark blue region to the right is ruled out by the requirement that perturbation theory 
remains valid up to the GUT scale. 

we calculate the spectrum of the Inert neutralinos, Inert charginos and Higgs bosons 
as well as their couplings, the branching ratios of the decays of the lightest CP-even 
Higgs state and the dark matter relic density. In order to calculate the dark matter 
relic density we use numerical methods. In particular, MicrOMEGAs 2.2 is used to 
numerically compute the present day density of dark matter. This includes the relevant 
(co-)annihilation channel cross sections and the LSP freeze-out temperature. MicrOMEGAs 
achieves this by calculating all of the relevant tree-level Feynman diagrams using CalcHEP. 
The CalcHEP model files for the considered model are generated using LanHEP pS]. The 
MicrOMEGAs relic density calculation assumes standard cosmology in which the LSP was 
in equilibrium with the photon at some time in the past. 
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6.1 Benchmark scenarios 



In order to construct benchmark scenarios that are consistent with cosmological obser- 
vations and colhder constraints we restrict our considerations to low values of tan/3 < 2. 
Figs. 1 and 2 show that in principle the appropriate value of dark matter density can be 
obtained even when tan/3 > 2. At the same time larger values of tan /3 lead to masses 
of the lightest and second lightest Inert neutralinos that are too small, as discussed in 
Section 3. As a result larger couplings of the lightest Inert neutralinos to Z are required 
to reproduce the measured value of f^coM^^- On the other hand according to Figs. 1 and 
2, light Inert neutralinos with substantial couplings to Z-boson give a considerable contri- 
bution to its invisible width leading to a conflict with LEP measurements (see discussion 
in subsection 6.2). 

However, even for tan/3 < 2 the lightest inert neutralino states can get appreciable 
masses only if either all or at least one of the Inert chargino mass eigenstates are light, 
i.e. m ± ~ 100 — 200 GeV. As clarified in Section 3 and in [26], the masses of the lightest 

Al 

inert neutralino states decrease with increasing m, ± and it is therefore rather difficult to 

Xl, 2 

find benchmark scenarios consistent with cosmological observations for m^± > 200 GeV. 
At the same time we demonstrate (see benchmark point (ix) in Table 2) that one light 
Inert chargino mass eigenstate is enough to ensure that the lightest inert neutralino state 
gains a mass of the order of M2/2. 

To obtain the kind of Inert neutralino and chargino spectrum discussed above one has 
to assume that some or all of the couplings Aq/j are rather small, e.g. they are expected to 
be much smaller than j^fi and /a/j. On the other hand in order to get m^o ~ m^o ~ Mz/2 
the Yukawa couplings fap and /^/j need to be relatively close to their theoretical upper 
bounds which are caused by the requirement of the validity of perturbation theory up to 
the GUT scale. Since gauge coupling unification determines the RG fiow and low energy 
value of g'l the mass of the Z' gauge boson is set by the VEV of the singlet field s only. 
In our study we choose s = 2400 GeV so that the Z' mass is about 890 GeV. This value 
of the Z' boson mass is just above the present lower bound of 865 GeV set by CDF [HUI 
and allows satisfaction of stringent limits on the Z' mass and Z — Z' mixing that come 
from precision EW tests |51] . 

Since we restrict our analysis to low values of tan /3 < 2 the mass of the SM-like Higgs 
boson is very sensitive to the choice of the coupling A. Stringent LEP constraints require 
A(Mt) to be larger than the low energy value of g[ ~ 0.47. If we try to increase X{Mt) much 
further, then the theoretical upper bounds on fa/3 and fa/s become substantially stronger. 
As a consequence, it is rather difficult to find solutions with \m^o\ ~ I"^x2l ~ 
Therefore in our analysis we concentrate on values of X{Mt) < 0.6. In addition, we 
set stop scalar masses to be equal to mq = mu = Ms = 700 GeV and restrict our 
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consideration to the so-called maximal mixing scenario when the stop mixing parameter 
Xt = At — \s / {^/2 tan (3) is equal to Xt = \/6Ms. This choice of parameters limits 
the range of variations of the lightest CP-even Higgs mass. In the leading two-loop 
approximation the mass of the SM-like Higgs boson varies from 115GeV(A = g[) to 
136GeV(A = 0.6). From Tables 1-2 one can see that the large values of A > (7^ that we 
choose in our analysis result in the extremely hierarchical structure of the Higgs spectrum, 
as pointed out in Section 4 (see also [23]). In Tables 1-2 the masses of the heavy Higgs 
states are computed in the leading one-loop approximation. In the case of the lightest 
Higgs boson mass the leading two-loop corrections are taken into account. 

The set of the benchmark points that we specify demonstrates that one can get a rea- 
sonable dark matter density consistent with the recent observations if \my^o \ ~ |m^o | ~ Mz/2. 
Our benchmark scenarios also indicate that in this case the SM-like Higgs boson decays 
predominantly into the lightest inert neutralinos (xi and X2) while the total branching 
ratio into SM particles varies from 2% to 4%. 

The benchmark points (i), (ii), (iv), (v) and (viii) are motivated by a non-Abelian 
family symmetry A27 which describes well the observed hierarchy in the quark and lepton 
sectors. As was discussed in Section 3 these scenarios imply that all flavour diagonal 
Yukawa couplings Xaa, faa and faa are rather small. Due to the approximate global U{1) 
symmetry fl2T|) . that originates from the family symmetry A27, the spectrum of Inert 
neutralinos involves a set of pseudo-Dirac states. When the masses of the lightest and 
second lightest Inert neutralinos are close or they form a Dirac state then the decays of 
hi into XaXiB "will not be observed at the LHC. Thus these decay channels give rise to a 
large invisible branching ratio of the SM-like Higgs boson. 

In Tables 1-2 we presented a few benchmark scenarios (i), (ii), (iv)-(vi), (ix) with 
almost degenerate lightest and second lightest Inert neutralinos. In some of these bench- 
mark points both lightest Inert neutralinos are lighter than Mz/2. Thus the Z-boson 
can decay into x«X/3 so that the lightest and second lightest Inert neutralino states con- 
tribute to the invisible Z-boson width. In other benchmark scenarios both of the lightest 
Inert neutralinos have masses above Mz/2 and the decays Z — )• XaXp are kinematically 
forbidden. 

When the LSP and NLSP are close in mass, LSP-NLSP co-annihilations may be an 
important factor in determining the dark matter relic density. If this is the case then 
the LSP-NLSP mass splitting should be an important factor. Since annihilations of two 
like-neutralinos are p-wave suppressed, one should compare (3Rzn with Rzu when trying 
to determine how important co-annihilations are, where /3 is the relative speed of the 
incoming particles, approximately 1/6. It is useful to consider the following situations. 
With the LSP and NLSP almost degenerate and with equal self-annihilation cross-sections. 
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but a negligible co-annihilation cross-section, the relic density of dark matter would be 
twice what it would have have been if the NLSP had not been present. If, alternatively, the 
co-annihilation cross-section was equal to the self-annihilation cross-sections, the existence 
of this extra channel would lead to a lower relic density. In this case it would in fact be 
equal to the relic density calculated in the absence of the NLSP. In this way, in such 
a scenario where co-annihilations and self-annihilations are about as important as each 
other, the rehc density is largely independent of the LSP-NLSP mass sphtting. 

For the benchmark scenarios (i) and (ii) this latter situation is approximately the case 
and the LSP-NLSP mass splitting turns out not to be an important factor. The mass 
splitting is in fact small, about half a GeV, but if it were larger and the NLSPs were 
made to have frozen-out much earlier, the relic density would only be decreased slightly 
(by about a tenth). In benchmark scenario (iv), even though the LSP and NLSP are 
close in mass, co-annihilations are unimportant due to the small value of Rzu- In this 
case increasing the NLSP mass substantially while keeping everything else fixed would 
lead to an approximate halving of the predicted relic density, since the NLSPs would 
have decoupled much earlier than, rather than at the same time as, the LSPs. The 
only other benchmark scenario where the LSP and NLSP are close enough in mass for 
co-annihilations to be potentially important is scenario (ix). Here co-annihilation is in 
fact the dominant process and changing the LSP-NLSP mass splitting would have a large 
effect on the predicted relic density. In fact, in this scenario, if the NLSP were not present 
the predicted relic density would be within the measured range. 

If the mass difference between the second hghtest and the lightest Inert neutrahno is 
10 GeV or more, then some of the decay products of a X2 that originates from a SM-hke 
Higgs boson decay might be observed at the LHC. In our analysis we assume that all 
scalar particles, except for the lightest Higgs boson, are heavy and that the couplings of 
the Inert neutralino states to quarks, leptons and their superpartners are relatively small. 
As a result the second lightest Inert neutralino decays into the lightest one and a fermion- 
antifermion pair mainly via a virtual Z. In our numerical analysis we did not manage to 
find any benchmark scenario with |m^o| — |m^o| > 20 GeV leading to reasonable values 
of l^cDM^^- Hence we do not expect any observable jets at the LHC associated with the 
decay of a X2 produced through a Higgs decay. However, it might be possible to detect 
some lepton-antilepton pairs that come from the decays hi — )■ X2Xa- In particular, we 
hope that fi^fi^ pairs that come from the exotic decays of the lightest CP-even Higgs 
state mentioned above can be observed at the LHC. 

In Tables 1-2 benchmark scenarios (iii), (vii), (viii) can lead to these relatively energetic 
muon pairs in the final state of the SM-like Higgs decays. Since the Higgs branching ratios 
into SM particles are rather suppressed, the decays of the lightest CP-even Higgs state 
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into + X might play an essential role in Higgs searches. 

In addition to the exotic Higgs decays, the scenarios considered here imply that at 
least two of the Inert neutralino states that are predominantly the fermion components 
of the Inert Higgs doublet superfields and one of the Inert chargino states should have 
masses below 200 GeV. Since these states are almost Inert Higgsinos they couple rather 
strongly to W and Z-bosons. Thus at hadron colliders the corresponding Inert neutralino 
and chargino states can be produced in pairs via off-shell W and Z-bosons. Since they 
are light their production cross sections at the LHC are not negligibly small. After being 
produced Inert neutralino and chargino states sequentially decay into the LSP and pairs 
of leptons and quarks resulting in distinct signatures that can be discovered at the LHC 
in the near future. 

6.2 Neutralino and chargino collider limits 

The remarkable signatures discussed above raise serious concerns that they could have 
already been observed at the Tevatron and/or even earlier at LEP. For example, the light 
Inert neutralino and chargino states could be produced at the Tevatron [52]. Recently, 
the CDF and DO collaborations set a stringent lower bound on chargino masses using 
searches for SUSY with a trilepton final state (i.e. trilepton signal) [53J. These searches 
ruled out chargino masses below 164 GeV. However this lower bound on the chargino 
mass was obtained by assuming that the corresponding chargino and neutralino states 
decay predominantly into the LSP and a pair of leptons. In our case, however, the Inert 
neutralino and chargino states are expected to decay via virtual Z and W exchange, i.e. 
they decay predominantly into the LSP and a pair of quarks. As a consequence the lower 
limit on the mass of charginos that is set by the Tevatron is not directly applicable to 
the benchmark scenarios that we consider here. Instead in our study we use the 95% C.L. 
lower limit on the chargino mass of about 100 GeV that was set by LEP II [54] . 

In principle LEP experiments also set constraints on the masses and couplings of 
neutral particles that interact with the Z-boson. As mentioned above when the masses 
of Xi and X2 are below Mz/2 they are almost degenerate and thus the decays of Z into 
XaXi3 contribute to the invisible width of the Z-boson changing the effective number of 
neutrino species N^^^ . The contribution of xi and X2 {^N^^^) to N^^^ is given by 

Amjf = (5n + 25i2 + ^22 , (38) 
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where 



8al3 — R%al3 



2M| ^ ^ M| 



2M| 




(39) 



Xa. I I I \ .1 Xa I I X^ I 



All three terms in Eq. f p8|) contribute to A^^-^-^ only if 2|m^o| < Mz- In the case where 
only the Z-boson decays into XiXi ^-^e kinematically allowed the values of 612 and 622 
should be set to zero. If l^^^-o] + 1^^X2 1 ^ while 2|m^o| > Mz then only and 612 
need to be taken into account. 

In order to compare the measured value of N^, with the effective number of neutrino 
species in the EgSSM, i.e. A"^-^-^ = 3 + AN^-I^^ , it is convenient to define the variable 

JV^ff — N^^P 

D = ^— , (40) 

where N'^^^ = 2.984 and a'^^^ = 0.008 [55]. The value of D represents the deviation 
between the predicted and measured effective number of neutrinos contributing to the Z- 
boson invisible width. It is worth pointing out that in the SM D = 2. In the benchmark 
scenarios presented in Tables 1-2 the value of D is always less than 3. Figs. 1 and 2 
also demonstrate that there is a substantial part of the EgSSM parameter space where 
m^o^ < Mz/2 and D < 3. This indicates that the relatively light Inert neutralinos with 
masses below Mz/2 are not ruled out by different constraints on the effective number 
of neutrinos set by LEP experiments (see, for example [55]-[56]). Indeed, as argued 
in Section 3 the Yukawa couplings fa/s and fai3 can be chosen such that the Rzai3 are 
very small. The couplings of the lightest and second lightest Inert neutralinos to the 
Z-boson are relatively small anyway because of the Inert singlino admixture in these 
states. Nevertheless Figs. 1 and 2 show that the scenarios with light Inert neutralinos 
which have masses below Mz/2 and relatively small couplings to the Z-boson can lead 
to the appropriate dark matter density consistent with the recent observations. 

LEP has set limits on the cross section of e~^e~ — )■ X2X1 ixtXi) the case when 
X2 QQXi ixf ~^ QQ'Xi) predominantly [57]. Unfortunately, the bounds are not directly 
applicable for our study because OPAL limits were set for a relatively heavy X2 ixf) 
only (l^^ol > 60GeV). Nevertheless, these bounds demonstrate that it was difficult 
to observe light neutralinos with l^^^-oj < 100 GeV if their production cross section 
a{e~^e~ — )■ XaX/?) ^ 0.1 — 0.3pb"\ Since at LEP energies the cross sections of colourless 
particle production through s-channel 7/Z exchange are typically a few picobarns the 
lightest and second lightest Inert neutralino states in the EgSSM could escape detection 
at LEP if their couplings Rzai3 ^ 0.1 — 0.3. 
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6.3 Dark matter direct detection 



Another constraint on the couphngs of the hghtest Inert neutrahno comes from exper- 
iments for the direct detection of dark matter. Recently the CDMSII and XENONIOO 
collaborations have set upper limits on the weakly interacting massive particle (WIMP)- 
nucleon elastic-scattering spin-independent cross section [l5],[58]. The XENONIOO Col- 
laboration claims a limit on the spin-independent cross section of 3.4 x 10~^^ cm^ for 
a 55 GeV WIMP. This limit remains fairly constant for lower WIMP masses and does 
not increase above about 4 x lO"'^*^ cm^ even for the lowest LSP masses that are con- 
sistent with our thermal freeze-out scenario. Since in the EgSSM the couplings of the 
lightest Inert neutralino to quarks (leptons) and squarks (sleptons) are suppressed, the 
X?-nucleon elastic scattering, which is associated with the spin-independent cross section, 
is mediated mainly by the t-channel lightest Higgs boson exchange. Thus in the leading 
approximation the spin-independent part of x'J^nucleon cross section in the EgSSM takes 
the form [591 ED] 



(41) 



where 



m^f^^ = {N\m,qq\N) , /^^ = 1 - ^ /^^ . 

q=u,d,s 

Here for simplicity we assume that the lightest Higgs state has the same couplings as the 
Higgs boson in the SM and ignore all contributions induced by heavy Higgs and squark 
exchangqj. Due to the hierarchical structure of the particle spectrum and the approximate 
Z2 symmetry this approximation works very well. Using the experimental limits set on 
asi and Eqs. (HTi) one can obtain upper bounds on X{\^ [62] . 

In Tables 1-2 we specify the interval of variations of asi for each benchmark scenario. 
As one can see from Eq. fl4ip the value of asi depends rather strongly on the hadronic 
matrix elements, i.e. the coefficients f^g, that are related to the vr-nucleon a term and 
the spin content of the nucleon. The hadronic uncertainties in the elastic scattering cross 
section of dark matter particles on nucleons were considered in [5^ [55] . In particular, it 
was pointed out that f^^ could vary over a wide range. In Tables 1-2 the lower limit on 
asi corresponds to f^^ = while the upper limit implies that = 0.36 (see [60]). From 



^The presence of almost degenerate lightest and second lightest Inert neutralinos could result in the 
inelastic scattering of Xi on nuclei {A), i.e. Xi + ^ ~^ X2 + that could affect the direct detection of 
Xi at the experiment. However such processes may take place only if the mass splitting between x? and 
X2 is less than 100 KeV [EI]. Since in all of the benchmark scenarios considered here the corresponding 
mass splitting is substantially larger the inelastic scattering of Xi does not play any significant role. 



26 



Tables 1-2 and Eq. ( HT]) it also becomes clear that asi decreases substantially when irih^ 
grows. 

Since in all of the benchmark scenarios presented in Tables 1-2 the lightest Inert 
neutralino is relatively heavy (|^7^;^o| ~ M^/2), allowing for a small enough dark matter 
relic density, the coupling of Xi to the lightest CP-even Higgs state is always large giving 
rise to a Xi^nucleon spin-independent cross section which is of the order of or larger than 
the experimental upper bound. However it is worth keeping in mind that the obtained 
experimental limits on a si are not very robust ^64J- Moreover, CDMS II and XENONIOO 
quote 90% CL. upper bounds while the 95% confidence level bounds are larger by a 
factor of 1.3 . By the same token the 99% CL. and 99.9% CL. upper bounds, which are 
associated with 2.6 and 3.3 standard deviations, are expected to be 2 and 3 times larger 
than the 90% CL. bounds respectively. Following these estimates it is clear that the 
benchmark scenarios presented in Tables 1-2 cannot yet be ruled out by either XENONIOO 
or CDMS II. However in the near future the expected new analysis from XENONIOO may 
either confirm or refute our scenario. 



7. Summary and Conclusions 

In this paper we have considered novel decays of the SM-like Higgs boson which can occur 
within a particular dark matter motivated scenario of the Exceptional Supersymmetric 
Standard Model (EgSSM). This model implies that at high energies the GUT gauge 
group is broken to the SM gauge group together with an additional U{1)m gauge group 
under which right-handed neutrinos have zero charge. To ensure anomaly cancellation and 
gauge coupling unification, the low energy matter content of the EgSSM includes three 
27 representations of Eq and a pair of SU{2) doublets from an additional 27' and 27'. 
Thus the EgSSM involves extra exotic matter beyond that of the MSSM that includes 
two families of Inert Higgs doublet superfields and and two Inert SM singlet 
superfields Sa that carry U{1)n charges. The fermion components of these superfields 
form Inert neutralino and chargino states. 

To satisfy LEP constraints we restricted our consideration to scenarios with relatively 
heavy Inert chargino states, i.e. m^± > 100 GeV. In our analysis we also required the 
validity of perturbation theory up to the GUT scale which sets stringent constraints on 
the values of the Yukawa couplings at low energies. Using these restrictions we argued 
that the lightest and the second lightest Inert neutralinos (x? and X2) always light, 
viz. they typically have masses below 60 — 65 GeV. These neutralinos are mixtures of 
Inert Higgsinos and singlinos. In our model Xi tends to be the LSP and can play the role 
of dark matter, while X2 tends to be the NLSP. The masses of Xi and Xo, can be induced 
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even if only one family of the Inert Higgsinos couples to the two SM singlinos. The masses 
of Xi and X2 decrease with increasing tan /3 and Inert chargino masses. 

An important requirement of this paper is that the lightest Inert neutralino account for 
all or most of the observed dark matter relic density. This sets another stringent constraint 
on the masses and couplings of Xi- Indeed, because the lightest Inert neutrahno states are 
almost Inert singlinos, their couplings to the gauge bosons, Higgs states, quarks (squarks) 
and leptons (sleptons) are rather small resulting in a relatively small annihilation cross 
section of XiXi ~^ SM particles and the possibility of an unacceptably large dark matter 
density. In the limit when all non-SM states except the Inert neutralinos and charginos are 
heavy (> TeV) a reasonable density of dark matter can be obtained for |m^o J ~ Mz/2 
where the Inert LSPs annihilate mainly through Z in the s-channel [26]. If Xi annihilation 
proceeds through the Z-boson resonance, i.e. 2|m^o| ^ Mz, then an appropriate value of 
^CDM^^ can be achieved even for a relatively small coupling of Xi to Z. 

The above scenario naturally emerges when a A27 family symmetry is included in 
the EgSSM [33]. The family symmetry was not introduced for this purpose, instead it 
was introduced earlier to provide an explanation of the Z2 symmetry and to account for 
the quark and lepton masses and mixings, including tri-bimaximal neutrino mixing. It is 
therefore encouraging to find that the same symmetry leads to a spectrum of inert pseudo- 
Dirac neutralinos which allows for a successful dark matter relic abundance, and also 
predicts novel Higgs decays. The A27 family symmetry also implies two almost degenerate 
families of Z)-fermion states [33] and in addition may have interesting consequences for 
S-physics [65]. As discussed in subsection (3.2) this symmetry leads to a cancellation of 
different contributions to the off-diagonal couplings of the LSP and NLSP. In addition, due 
to the singlino component of the lightest Inert neutralino states, the diagonal couplings 
of Xi and X2 to the Z-boson can also be rather small. Therefore these states could have 
escaped detection at LEP. 

The main point we make in this paper is that, within the above dark matter motivated 
scenario, although the lightest and the second lightest Inert neutralinos might have very 
small couphngs to the Z-boson, their couphngs to the SM-like Higgs state hi are always 
large. Indeed, we argued that in the first approximation the couplings of Xi and X2 to the 
lightest CP-even Higgs boson are proportional to |m^o J/ v. Since |m^o J ~ Mz/2 these 
couplings are much larger than the corresponding 6-quark coupling. Thus the SM-like 
Higgs boson decays predominantly into the lightest inert neutralino states and has very 
small branching ratios (2% — 4%) for decays into SM particles. We have illustrated this, 
together with the other phenomenological aspects of the dark matter motivated scenario 
considered in this paper, by presenting a set of benchmark points in Tables 1-2. If the 
masses of the lightest and second lightest Inert neutralinos are very close then the decays 



28 



of hi into XaXi3 "will not be observed at the LHC giving rise to a large invisible branching 
ratio of the SM-like Higgs boson. When the mass difference between the second lightest 
and the lightest Inert neutralinos is larger than 10 GeV the invisible branching ratio 
remains dominant but some of the decay products of X2 might be observed at the LHC. 
In particular, there is a chance that fi~^fi~ pairs could be detected. Since the branching 
ratios of hi into SM particles are extremely suppressed, the decays of the SM-like Higgs 
boson into l^l^ + X could be important for Higgs searches. 

In conclusion, the EgSSM predicts three Higgs families plus three Higgs singlets, where 
one family develop VEVs, while the remaining two which do not are called Inert. This 
pattern of Higgs VEVs is due to a broken Z2 symmetry whose origin can be understood 
from a A27 family symmetry. The model can account for the dark matter relic abun- 
dance if the two lightest Inert neutralinos, identified as the LSP and NLSP, have masses 
close to half the Z mass, with a pseudo-Dirac structure as predicted by the A27 family 
symmetry. Within this scenario we find that the usual SM-like Higgs boson decays more 
that 95% of the time into either LSPs or NLSPs, with the latter case producing a final 
state containing two soft leptons with an invariant mass less than or about 10 GeV. 
We have illustrated this with a set of benchmark points satisfying phenomenological con- 
straints and the WMAP dark matter relic abundance. This scenario also predicts other 
light Inert chargino and neutralino states below 200 GeV, and large LSP direct detection 
cross-sections close to current limits and observable soon at XENONIOO. 
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Tabic 1: Benchmark scenarios for m/^ « 133 — 135 GeV. The branching ratios and decay widths 
of the Ughtest liiggs boson, tlic; masses of the Higgs states, Inert neutrahnos and charginos as 
well as the couplings of Xi a-nd X2 ^^"^ calculated for s ~ 2400 GeV, A = 0.6, Ax = 1600 GeV, 
niQ = mu = Ms = 700 GeV, Xt = V6Ms that correspond to nih^ ~ Mz' ~ 890 GeV. 
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Table 2: Benchmark scenarios for m/^ « 114 — 116 GcV. The branching ratios and decay widths 
of the lightest Higgs boson, the masses of the Higgs states, Inert neutralinos and charginos as well 
as the couplings of Xi and x\ are calculated for s = 2400 GeV, A = g\ = 0.468, Ax = 600 GeV, 
niQ = mu = Ms = 700 GeV, Xt = VQMs that correspond to nih^ ~ Mz' ~ 890 GeV. 
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